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CANCELLATIONS IN POWER SERIES OF SINE TYPE 


J. ARIAS DE REYNA 


Abstract. We present a method to study the behavior of a power series of 
type 


/(*):= E(-l)" 


r 2n+l 


n =0 


(2n + 1)! 

when x —> oo. 

We apply our method to study the function 

rt rx j ny j 

f(t) := / — / — / —{sinrr + sin(a; — y) — sin(a: — z) — sin(:r — y + z)}. 

Jo x Jo y Jo z 

We will derive various different representations of f(t) by means of which it 
will be shown that limt—>.+oo /(£) = 0, disproving a conjecture by Z. Silagadze, 
claiming that this limit equals — 7r 3 /12. 


1. Introduction. 


In Titchmarsh m Section 14.32] we find two suggestive equivalents to the Rie- 
mann Hypothesis: The RH is equivalent to F{x) = 0{x^ +E ), or alternatively to 
G(x) = 0(:x~^ +e ) where 


( 1 . 1 ) 


Fix) = E 


(-l) r 


(n — l)!C(2n)' 


G{x) = E 


(- l) n x n 
nl((2n + 1) 


The first equivalence is due to Riesz [8j, and the second to Hardy and Littlewood 

m- 

Since £(n) converges to 1 these series can be considered slight modifications of 
the exponential function. As in the case of the exponential the series are convergent 
everywhere but the small values they get for x large is the result of an amazing 
cancellation between large terms of different signs. This phenomenon happens also 
in the case of the sine or cosine series. In these simple cases the many algebraic 
properties of the corresponding sums yield the proof of the cancellations. But how 


can one treat a case as the series (1.1) above? 


A problem in MathOverflow leads us to consider the power series 


2n+l 


( 1 . 2 ) 


/w = E<-i>"(E 


fc =i 


ffk-i 

k 


j.2n~hl 


(2n+l)!(2n+l) 


where the H n = A are the harmonic numbers. The proposer Z. Silagadze 

asked for a proof that /(f) —► —7 t 3 /12 (obtained by using some arguments from 
physics). 
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We think that our solution presented here is interesting because it provides an 
example of how to treat this type of problems. It must be said that the method 
can be applied to the functions in (1.11, but as we will see in Section [b] it only gives 
another path to the connection between these series and the Riemann Hypothesis. 
Silagadze’s problem in MathOverflow TO] was to compute 

r°° rtr ( x /I'll r v rl? 

(1.3) / -— / — / — {sin a; + sin(a; — y) — sin(:r — z) — sin(a: — y + z)} 

Jo x Jo V Jo z 

for which he conjectured the value —7r 3 /12. 

This integral is not absolutely convergent. Being a multiple integral it is not 
clear in what sense he is asking to compute it. The most natural interpretation of 

(1.3) is to define for t > 0 

... dx f x dy f y dz r . . . . , 

(1.4) /(f) := / — / — / —\smi + sm(i- y) — sm(i- z) — sm(x — y + z) 1. 

Jo x Jo V Jo z 

and define the integral (1.3) as the lim^oo /(f). We will prove that this limit = 0. 

The function /(f) defined by the triple integral in (1.4) has many interesting 
properties. It extends to an entire function whose power series (1.2) is a slight 
modification of the series for sint. We are interested in its behavior for t —> +oo. 
This is a similar problem as the one (equivalent to the Riemann Hypothesis) intro¬ 
duced by M. Riesz cited above. We will solve the problem in this relatively simple 
case, essentially, by means of a representation of /(f) as a Fourier transform (4.5). 

The transformations we apply to the function /(f) conclude with the represen¬ 
tation as a Fourier transform in Proposition |4.4| But all our representations are 
needed to reach this final one which will solve Silagadze’s problem. 

Our formulas (2.3) and (4.51 may also be applied to the computation of the triple 
integral. The numerical computation of a triple integral is quite often difficult. By 
means of the power series or the Fourier representation it can be computed easily 
and above all more reliably for relatively small values of f. For large values of f 
we present explicitly an asymptotic expansion (Section [5]) that is very well suited 
to compute /(f) with high precision for f large. The first term of the asymptotic 
expansion 


f(t) = 


cos f log“ f 


O 


(¥) 


2 f 

shows that /(f) has zeros near the zeros of cost for t large. 

As an application we will also obtain the x-ray of fit). It shows that the zeros 
with small absolute value are real. On average there are two zeros on each interval 
of length 27t. The first zeros are separated by approximately 4, 2, 4, 2, ... The x-ray 
also shows that the general properties of fit) are very similar to those of sint. 

To finish the paper we give some details about the Riesz function analogous to 
fit), also explaining why we may not hope that an analogous analysis will solve 
the RH. As expected! 


2. The entire function /(t). 

Proposition 2.1. The function fit) defined fort > 0 by the absolutely convergent 
integral 
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extends to an entire function with power series expansion 


2n+l 


(2.3) 


/» = £(-!)”( E 


H k - i 


^2n+l 


k ) (2n + l)!(2n+ 1) 


where H n = X]fe=i \ n-th harmonic number. 

Proof. Applying the Mean Value Theorem we get the following bound of the abso¬ 
lute value of the integrand in the definition of /(f) 

1 1 sinx — sin(x — z) sin(x — y + z) — sin(x — y) ^ 2 


x y 


xy 


When integrating over z £ (0, y) we get something bounded by 2/x, and integrating 
this over y £ (0, x) we get something bounded by 2. The integral of this over 
x £ (0,t) yields something bounded by 2 1. This shows that the integral in (2.2) is 
absolutely convergent. 

Changing variables x = tu, y = tv, z = tw yields 

[ l duf u dvf v dw, , 

/(f) = / — / — / — {sm(tu)+sm(tu — tv)—sm(tu — tw) — sm(tu — tv+tw)} 

Jo u J o v J o w 


and renaming the variables 
,-1 


,, . dx f x dy f y dz , . . 1 

/(f) = / — / — —lsm(fx) + sin(tx — ty) — sm(fx — tz) — sm(tx — ty+tz)\. 

Jo x Jo y Jo z 

This representation shows that /(f) extends to an entire function with power series 
expansion 

™ j.2n+l 


m = 

n—0 v ’ 


where 


A n = [ ldX [ Xd V f ^{x 2 " +1 + ( X _ yfn +1 _ (x _ ^+1 _ (x _ y + z) 2n+ 1} 
Jo x Jo 2/ Jo ^ 

Notice that A 0 = 0. To compute the other A n we get successively 

£ J {x 2 ” +1 - (* - z) 2n+1 }- = - XjV 1 )* ^ Y x2U+1 ~ k ’ 


dy 


dz 


2n+l 


^ / {x 2 " +1 -(x-z) 2n+1 }- = - ^(-1) 


/o V Jo 


k =1 


2n + 1\ x 2n+1 

~W~ 


dx 


dy 


. .. 2n+l 

- (x= -—E(-i) 


“ / ^ / { 
/o ^ Jo 2/ Jo 


2n+l 


2n + 1 


fc=i 


2n + l\ 1 
k 

2n+l 


2n 


1 1 " 

— y 
+ 1 ^ 

fc=l 


fc 2 

Hk 

k 


where the last step is proved in Lemma [2~~i| 
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For the other part of A n we compute 


j\(x - yf "« - (x - „ + x) 2 "«}f = - g ( 2 " + y j(x - 

f ^ [ V {( x ~ y) 2n+1 - (x - y + x) 2n+1 }— 

Jo y Jo z 


2n+l 


= -£ 


1 /2n + 1 


fe=i 


) L yk ^ x ~ ^ 


2ra+l-fc dy 

y ' 


In the last integral we change variables y = xu yielding 

, dit 


x 2n+1 (2n + 1\ /- 1 


/c=l 




2n+l—/c 1 


y; 1 ® 2 ”-*' 1 (2n + 1)! (k — l)!(2n + 1 — fc)! _ y* 1 a; 2n+1 

k k\(2n +1 — k)\ (2n+l)! 


fe=i 


fc=i 


k 2 


so that 


1 v f ? / V - ») 2 " +1 -<*-»+ 2 > 2 ” +, }f =-^1;^ 

^ fe=i 


^ 2n+l . 


do a-’ do y JO 
Collecting all our partial results we obtain 

1 2Ti+l „ 2n+l 1 271+1 „ 

^ _ t yi rife 1 yi 1 _ 1 yi Hfc-1 

n 2n + 1 2n + 1 fc 2 2n + 1 

k— 1 fc=l 


□ 


Lemma 2.4. For any natural number n we have 

ic=i v 7 fc=i 


Proof. To prove a n = b n it is sufficient to prove a\ = &i and a n — a„_i = b n — & n _i- 
In our case the equality for n = 1 is checked easily. So, we have to show 


£(-i ) fc+1 

k =1 


1_ 

k 2 \ k ) k 2 n 

k =l v 7 


Since (^) = (^_J) + ("fe 1 ), the equality is equivalent to 


(~ir +1 , v f _nHi 

n 2 ' U-l/fc 2 n 

fc=l V 7 

or, multiplying by n 

m X >‘>‘ +1 (1)1 = ^ 

fc = l X 7 

This is formula 0.155.4 in Gradshteyn and Ryzhik 0 P- 4]. □ 
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3. Plot of the function f(t). 


The power series expansion of fit) is analogous to the one for sinf. These 
power series are not well suited for computation, because of the violent cancellation 
between large terms. Nevertheless, one may use it using high precision in the 
computation of the terms to get approximate values. Assuming t > 0, each term of 
the power series (2.3) is in absolute value less than t 2n+1 /(2n + 1)!. Therefore the 
error committed by summing only the first N > 3t terms of the series (2.3) is less 
than 2~ 2N . It follows that to compute f(t) with error less than 2e we need only 
compute the sum of the first N terms of the series with error less than e, taking N 
large enough so that 

log(l/e) 


N > 


2 log 2 


N > 3t. 


All the terms of the series are less than t 2n+1 /{2n + 1)! < ( 2 n + 1 ) 2n+1 < e*. So, we 
must compute each term with an error less than e/N, for which it will suffice to 
compute each term working with a precision 


„ t + log(N/e) , . , . 

P := ----- decimal digits. 

log 10 

Of course this will be difficult for t very large, but today we may easily compute 
with thousands of digits of precision. 


0.5 

0.0 


-0.5 


- 1.0 


-1.5 


i A MA/WV^^ 

V 


We may say that the computation of the power series is much easier than the 
computation of the multiple integral (2.2) or any other multiple integral giving f{t) 
considered in this paper. This is true even when t is small. 


4. Computation of the limit of f(t). 
First we prove the following integral representation of f(t): 

Proposition 4.1. For t > 0 we have 


(4.2) m = 


1 r l 


du 


x sm u 


sin(ztx) 


sin (uxy) 


/o J o W(1 - ar)(l - xy) u(l-x)(l-y) u(l - y)(l - xy) 
Proof. Starting from the power series expansion it is easy to get 

£ 271+1 

^ k J (2n + l)!(2n+ 1) 


dx dy. 


m = £(-!)" E j £ 7 


71 = 1 


2ti+1 


k -1 


= du dx dyJ2(- l) n ^ x k 




,2 n 


71=1 


k =1 


3 =1 


(2n + 1)! ’ 
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If C is a circle with radius r > t, we may express the factorial by means of the 

Residues Theorem 

(4.3) 


/(f) = / du 


dx 


1 « 2 oo 2n+l k— 1 

„„ „„ 

•'ll jo n=1 k=1 j—r 

(Since |«/z| < 1 all these series converge absolutely and the order of the integrals 
and summations may be interchanged). This can be written as 


/(f) = 


du / dx 


dy- 


1 


7 4 3 o 24 224 224 

e z u x yz* — u z xz* — u z x z z* — u z x z yz* 


2ni J c z 2 (u 2 + z 2 ){u 2 x 2 + z 2 ){u 2 x 2 y 2 + 2 2 ) 

By the Residue Theorem we obtain 


dz. 


1 

27ri 


^4 3 2 24 224 224 

e z u x yz* — u z xz* — u z x z z* — u z x z yz* 


l c z 2 (u 2 + z 2 )(u 2 x 2 + z 2 )(u 2 x 2 y 2 + z 2 ) 


dz 


x sin u 


sin(m’) 


+ 


sin(ua:y) 


u(l-x)(l-xy) u(l-x){l-y) u(l - y)(l - xy) 


establishing our claim (4.2). Observe that since our integral representation (4.2) 
appears after computing one of the integrals in the absolutely convergent integral 
(4.3), the integrals in our new representation are also absolutely convergent. □ 


Notice that in the representation (4.2) we may interchange the integral in u with 
the integral in (x, y ), but then it is not easy to justify the interchange of the limit 
in t with the integrals, because the integrand is not dominated by an integrable 
function. If we would proceed formally in this way we would easily get that the 
limit is 0, but, as said before, this is not allowed. Therefore we follow another path. 


Proposition 4.4. For any complex t we have 


dx. 


n— 1 


Proof. We take t > 0 in the representation (|4.2[) , and change the order of integration 


/(*) = 


dx dy 


/o 


x sm u 


sin(iiir) 


sin (uxy) 


i(l-x)(l-xy) u(l-x)(l-y) u(l - y)(l - xy) 


du. 


Now subdivide the inner integral in three and change variables appropriately 


/(f) = 


dx dy 


o Jo 


xsmu 


0 u(l-x)(l-xy) 
f tx sinrt 


du 


r t.xy 


I o u{l-x)(l-y) 


du - 


sm u 


u (! - y) 0 - - x v) 


du 


We use Iverson’s notation [5], so that, for any proposition P, the symbol [P] is 1 if 
P is true and 0 if it is false. In this way we may write 


/(f) = 


dx dy 


o 


x[u < t] sin u ^ \ u < f*] si nu 


[u < txy] sin u 


i(l - x){l - xy) 


(l-x){l-y) u(l - y)(l - xy). 


du. 
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To simplify the notation we put a := u/t , so that always 0 < a < 1. Interchanging 
the integrals we obtain 


/(*) = 

r t 


sin u 


du 


/o 


x 


du—~ 


\a < x\ 


[a < xy] 


lo Jo \{l-x)(l-xy) (l-x)(l - y) (1 - y){l - xy) 

and will compute the inner double integral in x and y. This is 

rl rl ' ~ r ~ ^ '" 1 r a < xy] 


dx dy 


J{a) : = L L ( 


x 


, fa < x] 

du - „ L ,,, J —- + 


^ dxdy. 


I o Jo \(l~x)(l-xy) (l-x){l-y) {l-y){l-xy). 

We subdivide the square (0, l) 2 into three disjoint sets 

Si := {(x, y) : 0 < x < a, 0 < y < 1}, S 2 := {(x, y) : a < x < 1, 0 < xy < a}, 

S 3 ■= {(%,y) ■ a < xy}. 

In Si the integrand is = ^ 1 _ x ^ 1 _ xy ^ j , and in S 2 it is 

x 1 1 

(l-x)(l-xy) (l-x)(l-y) (1 2/)(l xy) ‘ 

In S 3 the integrand is equal to 

x 1 1 


(l-x)(l-xy) (1 — x)(l — y) (l-y)(l-xy) 


= 0 


so that 


J(a) = r dy f 
Jo Jo 


nd/x 


0 Jo (l-x)(l-xy) 


dx dx 


1 


0 (1 - Z/)(l - xy) 

1 ^2 

= -log 2 (l-a)- — + £ 


dy 

(1 -ay 


6 z —' n* 

n= 1 


This yields 


, f l ( 1. 2 / (1 — u/t) n — 1\ sina 

m = l (j >°g (i - »/*) + E -—U —) at du ■ 

A change of variables u = tx yields (4.5) for t > 0, and it is clear that the two sides 
are entire functions. □ 

Proposition 4.6. We have 

1,1 1 °° (1 — x) n — In sin(ta) 


/■ /I 

lim /(f) = lim / ( - log 2 (l - x) + ^ ■ 

>-oo t —>-+oo / n \ Z z ' 

n—1 


dx = 0. 


Proof. This is an example of the Riemann-Lebesgue Lemma, once it has been shown 
that the function 


g(x) := + 

~ n—1 

is in jC 1 [0,1]. This is a simple exercise. 


(l-x) n ~l 


□ 
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5. Asymptotic Expansion. 


The representation (4.4) as a Fourier integral yields by known methods an as¬ 
ymptotic expansion for /(f). First notice that the dilogarithm function Li 2 (;r) can 
be defined in the cut plane C\ [1, oo) by 


(5.1) 


Li 2 (z) = - [ log(l - u) 
J o 


du 

u 


where the path of integration is the segment joining 0 and z. For z in the plane 
with two cuts along the real axis (—oo, 0] and [1, +oo) the dilogarithm satisfies the 
Euler functional equation (cf. Lewin [BJ (1.12), p. 5]) 


(5.2) 


Li 2 (^) + Li 2 (l — z) = — - logzlog(l - z). 
6 


Here and in the sequel we will denote by log w the principal branch of the logarithm. 

Definition 5.3. Let SI be the complex plane with the two cuts (—oo, 0] and [1, +oo). 
We define the function g[z ) holomorphic in S7 by 


(5.4) g(z) := (- log 2 (l - z) - y + Li 2 (l - zfj - 

= log 2 (l -z)- Li 2 (z) - log z log(l - zfj^. 

It is easier to give explicitly the asymptotic expansion of the Fourier transform 

(5.5) 


J(t) := f g(x)e ixt dx. 
Jo 


Since fit) = Im J(t) its asymptotic expansion can be obtained easily from the one 
of Jit). 

Proposition 5.6. For any t > 0 we have f{t) = Im J(t) where 

pOO pOO 

(5.7) J(t)=i g[iy)e~ ty dy — ie lt / g{\ + iy)e~ ty dy. 

Jo Jo 


Proof. By the definition of g{z) and (4.4) we have 


f{t) = Im f gix)e lxt dx. 
Jo 


We apply Cauchy’s Theorem to a rectangle with vertices at 0, 1, 1 + iR, iR with 
R > 0 and let R —t +oo. In this way the integral in (0,1) can be converted into 
the two infinite integrals in the Proposition. The bounds are easy. □ 


Since the two integrals in Proposition 5.6 are Laplace integrals we may apply 
Watson’s Lemma [3, 17.03, p. 501] to get their asymptotic expansions. In this 
case we have the additional difficulty of logarithmic singularities at the extremes of 
the integrals at 0 and 1. We follow the path in Lyness [7]. To get the asymptotic 
expansion we need the behavior of g{z) near z = 0 and z = 1. 
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Lemma 5.8. For \z\ < 1 in Q we have 
' H n+ 1 2 


(5.9) g{z) = E 


n—0 


n + 1 (n + iy 


n—0 


oo oo 

:=Y A nZ n +(Y R n Zn ) 1 o S 2 - 


n—0 


n—0 


For z £ Q, with |1 — z\ < 1 we have 

OO 1 OO 

(5.10) g(z) = Y *'(" + 0(1 - z) n + - (Y0 - *)") l«g 2 (l - z) 

n—0 

oo oo 

“ E - z) n + (E DnO - z) n ) log 2 (l - z) 


n—0 


n—0 


n—0 


where \J r ( 2 ) = V(z)/T(z) is the digamma function. 


Proof. Near z = 0 we have g(z) = gi{z) + < 72 (z) log 2 where by (5.4) we have 

. l 0 g(l- 2 ) 


9i(z) = log 2 (l - z) - Li 2 ( 2 )) g 2 (z) 


< 71 ( 2 ) and < 72 ( 2 ) are holomorphic at 2 = 0. Expanding in power series at 2 = 0 we 
get fl5Jl ). 

Near 2 = 1 we have g(z) = 53 ( 2 ) + < 74 ( 2 ) log 2 (l — 2 ) with < 73 ( 2 ) and 94 ( 2 ) holo¬ 
morphic at 2 = 1. By (5.4) 

7 r 2 Li 2 (l — 2 ) . . 1 

S3<2) = -fc +- 1 -' 9l(s) = 2 ? 

The functions ( 73 ( 2 ) and ( 74 ( 2 ) are holomorphic at 2 = 1 and expanding them in 
power series we obtain (|5.10). We notice the equalities. 


(5.11) 


V(n+1) = H n -j, ^'(^ +!) = tt - E fc 2 ‘ 


k —1 


□ 


Proposition 5.12. The following asymptotic expansion is valid for t -A +00 

OO | OO | 

(5.13) j(f)~E* n+1 ^^ + E* n+1Sra (^( n + 1 ) _logi + y) 


n—0 


n—0 


t n+l 


~ ielt Y^ nCn ^+r~ ielt E( _ 0"- D n{^'( n + 1 )+(7-^+logt+y) | 

n—0 n—0 

where the coefficients A n , B n , C n and D n are defined in Lemma\5. 


i! 


t n+l 


Proof. By the results in Lyness [7] the asymptotic expansion of J(t) is obtained by 
integrating term by term the expression (5.7) after substituting the power series at 
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the extremes at 0 and 1, respectively. Explicitly 


m ~ x] * n+1 A » 

n =0 
oo 

-ie"^-*)"^ j 


y n e~ tv 


00 /*oo . • 

dy+^i n+1 B n y n (log y + — J e~ ty dy 


n—0 
co 


n—0 


y n e~ t ydy-id it Y J {-i) n Dn 

n —0 


y- y) e 


Computing the integrals we obtain (5.13). 


□ 


Corollary 5.14. The first order terms of the asymptotic expansion of f{t) are 
(5.15) 


cos t log" t | y sin t 
~2 ~ + V 2 


— 7 cos i — 1 


)¥+( 


57r 2 cos t y 2 cos t 77r sin t \ 1 

—4 - 2 1 2 — 1 -' |, h 


where 7 is Euler’s constant. 


6. Riesz type functions 


We have proved that the function f(t) given by the power series (1.2) tends to 
0 by writing it as a Fourier transform of an L 1 function. As an example of what 
happens in the case of a Riesz-type function we just consider Riesz’s F(x). Riesz 
proves [Sj 


( 6 . 1 ) 


F(l > = £ tTz 


n+l^n 


n =1 


(n- l)!C( 2 n) 


= *£ 


p{n) 


The Riemann Hypothesis is equivalent to F(x) = 0(xi +E ) for any e > 0. Riesz 
shows that F(x) is not 0(x a ) for any a < 1/4. Therefore, it is not true that F(x) 
converges to 0 when x —> 00 . Riesz proves that F{x)x~ x / 2 converges to 0. With 
some work we may reprove this. Noticing that 


(6.2) Vt e~* 

we can prove 


1 f +oa cos(§ arctamr) 

7*1 77*W‘ ™ {xt)dx ' 


t > 0 


(6.3) 





ncos(§ arctan(rc 2 i)) \ 
(LkW (i + „.(2)3/. J 


cos(a :t) dt. 


We apply summation by parts and use M(x) := J/ n <„ yfn) 
show that 


(6.4) 


OO 

g(.t) ~ £ 

n—1 


ncos(| arctan (n 2 t)) 

(1 + n 4 t 2 ) 3 / 4 


0( xe 


to 


is in L 1 (0, 00 ). This implies that F(x)x~ 1 ^ 2 = o(l). 

We do not give the details of the proofs, but they are clearly more complicated 
than those given by Riesz. 

We can also obtain in this way that M(x) = 0{x a ) with \ < a < 1 would imply 
that £(s) does not vanish on Re s > a. Again quite a difficult way to get this simple 
result. 
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This represents f(z) on the square (—2,30) x (—16,16). In red lines where f(z) 
takes real values, in blue lines where f(z) is purely imaginary. At t = 0 we observe 
a triple zero. Table 1 contains some more zeros of f(t). 
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Table 1. First zeros of f ( t ). 


1 

7.67705050991057 

16 

53.7503680917393 

31 

101.945883830959 

2 

8.93841793140833 

17 

57.9542719238923 

32 

104.179794651601 

3 

13.9549229413771 

18 

60.0643806701981 

33 

108.230057788273 

4 

15.5679745247884 

19 

64.2391371546902 

34 

110.475488483408 

5 

20.2405336746206 

20 

66.3740436665322 

35 

114.514173951783 

6 

22.0185470304095 

21 

70.5238571367734 

36 

116.770227743121 

7 

26.5267465905312 

22 

72.6803262355070 

37 

120.798238189933 

8 

28.4064753351502 

23 

76.8084524955294 

38 

123.064133569705 

9 

32.8127796622118 

24 

78.9839168600548 

39 

127.082255551763 

10 

34.7636137360244 

25 

83.0929401032590 

40 

129.357306301363 

11 

39.0985258618221 

26 

85.2853203526836 

41 

133.366230408727 

12 

41.1028595266668 

27 

89.3773338407949 

42 

135.649829884233 

13 

45.3839993880412 

28 

91.5849167023848 

43 

139.650166567860 

14 

47.4305876227269 

29 

95.6616452128959 

44 

141.941775186108 

15 

51.6692361998946 

30 

97.8829983472283 

45 

145.934067362779 


[5] D. E. Knuth, Two notes on notation , Amer. Math. Monthly 99 (1992) 403-422. http:// 
arxiv.org/abs/math/9205211 

[6] L. Lewin, Poly logarithms and Associated Functions , North Holland, New York, 1981. 

[7] J. N. Lyness, Adjusted Forms of the Fourier Coefficient Asymptotic Expansion and Applica¬ 
tions in Numerical Quadrature , Math. Comp., 25 (1971) 87-104. 

[8] M. Riesz, Sur Vhypothse de Riemann , Acta Math. 40 (1916) 185-190. 

[9] Z. K. Silagadze, Multiple Integral (American Mathematical Monthly problem 11621 and its 
generalizations), http://mathoverflow.net/questions/125237 March 22, 2015 

[10] Z. K. Silagadze, Interesting triple integral , http://mathoverflow.net/questions/200880 
March 24, 2015. 

[11] R. Stong, A quadruple integral , Amer. Math. Monthly 120 (2013) 667-669. 

[12] E. C. Titchmarsh, The theory of the Riemann zeta-function, 2nd ed. with a preface by D. R. 
Heath-Brown, Oxford University Press, New York, 1986. 

Facultad de Matematicas, Univ. de Sevilla, Apdo. 1160, 41080-Sevilla, Spain 
E-mail address: arias@us.es 










